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Abstract 

Heterosis is associated with differential gene expression between hybrids and their parental lines, and the genes involved in 
cell proliferation played important roles. AtARF2 is a general cell proliferation repressor in Arabidopsis. In our previous study, 
two homologues (ZmARFlO and ZmARF25) of AtARF2 were identified in maize, but their relationship with heterosis was not 
elucidated. Here, the expression patterns of ZmARFlO and ZmARF25 in seedling leaves of maize hybrids and their parental 
lines were analyzed. The results of qRT-PCR exhibited that ZmARF25 was down-regulated in leaf basal region of hybrids. 
Moreover, overexpression of ZmARF25 led to reduced organ size in Arabidopsis, which was mainly due to the decrease in cell 
number, not cell size. In addition, the cell proliferation related genes AtANT, AtGIFI and AtCRF5 were down-regulated in 
35S:\ZmARF25 transgenic lines. Collectively, we proposed that the down-regulation of ZmARF25 in maize hybrid may 
accelerate cell proliferation and promote leaf development, which, in turn, contributes to the observed leaf size heterosis in 
maize. 
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Introduction 

Heterosis is defined as advantageous quantitative and qualita- 
tive traits of offspring over their parents, and the utilization of 
heterosis principles has been a major practice for increasing 
productivity of plants and animals [1]. In spite of its importance 
and successful application, the underlying molecular mechanism 
for heterosis is still needs to be elucidated [2-5] . Studies on plant 
form evolution suggest that phenotype evolution often proceeds 
through changes in the spatial and temporal patterns of gene 
expression [6] . Theoretically, all the genes in hybrid are inherited 
from its two parental lines, but hybrid performance or phenotype 
can be quite different from its parents, thereby, demonstrating 
heterosis [7,8]. Therefore, the alteration of gene expression in 
hybrids may be responsible for hybrid vigor [4,8]. Up to date, 
genome wide gene expression profiles of different organs or tissues 
between hybrids and their parents in plants have been analyzed 
[9-15], and all possible modes of gene action were observed, 
including additivity, high- and low-parent dominance, under- 
dominance, and overdominance [9-15]. To gain further insight 
into the possible roles of differentially expressed genes between 
hybrid and its parental lines in plant development and heterosis, 
the functions of some differentially expressed genes have been 



characterized, and some of which may contribute to the observed 
heterosis. For example, overexpression of wheat hybrid up- 
regulated gene (ADP-ribosylation factor, TaARF) in Arabidopsis 
displayed increased growth rates during the seedling and adult 
stages [16]. Overexpression of Larix hybrid up-regulated gene 
LaAP2Ll, a member of the AP2/EREBP transcription factor 
family, led to markedly enlarged organs and heterosis-like traits in 
Arabidopsis [17]. Expression of an up-regulated protein ZmACT2 
in dry embryos of maize hybrid in Arabidopsis partially comple- 
mented the low germination phenotype in the Atact7 mutant [18]. 

Heterosis affects nearly any trait in nearly every hybrid, hence 
one of the major challenges for expression analysis is deciding 
which tissue to be analyzed [19,20]. It was proposed that seedling 
leaves can be used as a model system to investigate the molecular 
basis of maize heterosis, because it is commonly observed that 
hybrids produce larger leaf areas than their corresponding 
parental lines [20-24]. At the cellular level, a crucial effect of 
leaf size heterosis is manifested by primarily increases in cell 
number, not cell size [25-27]. Recent study provided molecular 
evidence for the relationship between cell proliferation and 
heterotic maize growth. For example, the expression level of 
maize CNR gene, ortholog of fruit weight 2.2 (fw2.2) in tomato, was 
founded to be negatively correlated with growth activity and 
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Figure 1. The third leaf areas and expression patterns of ZmARFlO and ZmARF25 between maize hybrids and their parental lines. (A) 

The third leaf area and mid-parent heterosis at the 8 DAG. H represents mid-parent heterosis (MPH). **refers to significant at p<0.01, data are 
represented as the mean ± SE (n = 3). (B-D) Expression patterns of ZmARFlO and ZmARF25 in different leaf regions of maize hybrids and their 
parental lines at the 5 DAG. B, C and D indicated as basal region, tip region and basal region, respectively, fi-actin was used as an endogenous control. 
Data are represented as the mean ± SE (n = 3), *refers to significant at p<0.05, **refers to significant at p<0.01. 
doi:1 0.1 371 /journal.pone.0094830.g001 



hybrid vigor [28]. Collectively, it reveals that the study of cell 
number control is eminently important to understanding the 
underlying mechanism for leaf size heterosis. 

In Arabidopsis, auxin signaling plays a vital role in leaf growth 
and development processes and some key genes, including 
AtARF2, AtARF7 and AtARF19, have been identified [29-31]. In 
our previous study, 31 maize genes that encode ARFs (Auxin 
Response Factors) were identified in maize genome, and 
phylogenetic analysis demonstrated that two of them, ZjnARFlO 
and ZjnARF25, were homologous to Arabidopsis growth repressor 
AtARF2 [32]. In this study, the expression patterns of %mARFl 0 
and ZjnARF25 in seedling leaves of maize hybrids and their 
parental lines were investigated and the function of hybrid down- 
regulated gene ZjnARF25 was characterized. Our data showed that 
ectopic expression of ZjnARF25 in Arabidopsis led to reduced organ 
size, indicating that the alteration of auxin signaling pathway 
might play an important role in maize leaf size heterosis. 



Materials and Methods 

Plant Materials and Growth Conditions 

Maize inbred lines Mo 17, B73, Zong3, 87-1 and their hybrids 
Mol7/B73, B73/Mol7, Zong3/87-l, 87-l/Zong3, B73/Zong3, 
Zong3/B73, B73/87-1 and 87-1/B73 were used for this study. 
Maize seeds were imbibed at 28°C in the dark, and the seeds after 
germination were selected with the same length of coleoptile 
(2.0 cm) to grow in hydroponic culture on a 16 h light/8 h dark 
cycle, at a temperature of 29°C and 25°C for the light and dark 
cycle, respectively. In maize leaf, the mature zone is composed of 
fully differentiated cells and cell division is restricted to the growing 
zone composed of immature cells [19,33]. Therefore, thirty 
millimeters from the base of the third leaves and twenty 
millimeters from the tip of the same leaf of a seedling 5 days 
after germination (DAG) were dissected and harvested for RNA 
extractions. 
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Figure 2. Amino acid sequences alignment and graphical representation for domain organization of ZmARF25 proteins. Identical 
amino acids are indicated by asterisks, amino acid position numbers are indicated on the right. B3 domain is marked with green box, ARF subdomain 
is marked with blue box, MR is marked with red box, and subdomain III and IV are marked with purple box and orange box, respectively. 
doi:1 0.1 371 /journal.pone.0094830.g002 



Arabidopsis (Col-0) seeds were sterilized by a 60 s 70% ethanol 
treatment followed by 1 % NaCIO within 1 0 min, four washes in 
distilled water, and chilled at 4°C for 3 days. Seeds were sown on 
half-strength Murashige and Skoog (MS) plates with 3% sucrose 
and 0.8% agar in a culture room at 22°C under a 16/8 light/ dark 
cycle. After germination, 10-day-old Arabidopsis seedlings were 
transplanted and grown at a density of four plants per pot 
containing a mixture of soil and vermiculite (2:1) at 22°C under a 
16/8 light/dark cycle with 70% relative humidity. 

Leaf Area Measurements 

Leaves were scanned using ImageScanner (GE Healthcare, 
USA), and measured with a public domain image analysis software 
(ImageJ version 1.32; http://rsb.info.nih.gov/ij/). For maize leaf 
area measurements, the third leaves were harvested in the 
morning (10:00-10:30 am) and experiments were done in three 
biological replicates and 1 0 plants of each genotype were analyzed 
for each replicate. The f-test was used to determine the significance 
of the differences in leaf areas between hybrid and the 
corresponding mid-parental value. The mid-parent heterosis 
(MPH) was calculated using the following formula: MPH = (Fi — 
mid-parental value) /mid-parental value in %, where F] is the 
average value of hybrid, and the mid-parental value is the average 
value of the two parents. 

Ten 30-day-old Arabidopsis plants of each transgenic line and 
wild-type were selected randomly for the fully expanded fifth leaf 
area measurements. The it-test was used to determine the 
significance of the differences in leaf area between transgenic 
plants and wild-type plants. 



Histological Analysis 

After the area measurements, the fully expanded fifth leaves of 
the ten 30-day-old Arabidopsis plants of each transgenic line and 
wild-type were used for epidermal cell size measurements. Leaves 
were first boiled in water for 10 min to kill the cells, and then 
immersed in 95% ethyl alcohol, after decolorized, they were 
transferred to hot (96°C) 85% lactic acid for 8 min. Finally, the 
leaves were transferred to room temperature lactic acid and 
examined using a Nikon Stereo Microscope. The number of 
epidermal cells was determined in an area of 0.58 mm in the 
middle region of the fifth leaf blade. The average cell size was 
measured by dividing 0.58 mm by the number of epidermal. The 
total number of cells per leaf was calculated by dividing the leaf 
area by the average cell size. Statistical calculations were 
performed with Microsoft EXCEL. 

Isolation of DNA and Total RNA 

Total DNA was extracted by the cetyltrimethylammonium 
bromide (CTAB) method [34] with minor modifications. Total 
RNA was extracted using a standard Trizol RNA isolation 
protocol (Invitrogen, USA), and then treated with RNase-free 
DNasel (Promega, USA) to remove genomic DNA. The amount 
and quality of the total RNA were verified by electrophoresis on a 
1 % agarose gel. The concentration of RNA was measured using a 
NanoDrop spectrophotometer (ND- 1000; NanoDrop Technolo- 
gies). First-strand cDNA synthesis was performed using 2 u.g of 
DNase digested total RNA with oligodT15 primer according to 
the manufacturer's protocol for RT-PCR first-strand synthesis 
(TaKaRa, Japan). 
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Figure 3. Overexpression of ZmARF25 in Arabidopsis. (A) The 

genomic PCR results for four homozygous T 3 transgenic plants 
(3SS::Zrr\ARF25 #1-4) using Bar gene specific primer. PCR products 
were determined visually by running on 1% agarose gel and stained by 
Ethidium bromide. Col represents wild-type and #1-4 represent 
35S::ZmARF25 transgenic lines, which is the same in Figure 3B; "+" 
represents the recombinant plasmid. (B) Relative expression levels of 
ZmARF25 gene in four homozygous T 3 transgenic lines. Total RNA was 
isolated from the 30-day-old aerial of Arabidopsis. Data are represented 
as the mean ± SE (n = 3). ACTIN2 was used as an endogenous control. 
doi:1 0.1 371 /journal.pone.0094830.g003 

Quantitative RT-PCR Analysis 

Quantitative RT-PCR (qRT-PCR) using SYBR Green PCR 
master mix (TaKaRa, Japan) was performed using CFX96 Real- 
Time PCR Detection System (Bio-Rad Laboratories, Inc., USA). 
Specific primers for qRT-PCR analysis were listed in Table SI. 
The PCR efficiency of all primer pairs was determined from 
standard curve experiments. PCR conditions consisted of an initial 
step at 95 °C for 3 min followed by 40 cycles of 95 °C for 15 s, 
followed by 60°C for 15 s, and then 72°C for 30 s. For 
amplification product specificity, a melting curve was generated 
at the end of each run, and verified by agarose gel electrophoresis 
of PCR products. All reactions were run in triplicate, and each 
reaction was completed with a melting curve and amplification 
curve analysis to confirm the amplification specificity. Ct values 
were determined using the CFX96 software with default settings. 
Differences between the Ct values of target gene and Actin were 
calculated as ACt = Ct target gene — Ct^ c(1 „, and the relative expres- 
sion levels of target genes were determined as 2 _ACt . For each 
sample, PCR was performed on three biological replicates. The 
average values of 2 _ACt were used to determine difference in gene 
expression. Actin was amplified as an endogenous control. 
Statistical analysis of the difference in the relative expression 
levels of target genes was performed by using Hest. 

Gene Cloning and Sequence Analysis 

PCR was performed using the ZmARF25-F/R primer pair 
from the first-strand cDNA from the B73 seedling leaves. PCR 
conditions consisted of an initial step at 94°C for 5 min followed 



by 35 cycles of 94°C for 25 s, 58°C for 25 s, and 72°C for 2 m 
45 s; and then a final step at 72°C for 10 min. The PCR products 
were confirmed by visualization of the samples by electrophoresis 
on 1 % agarose gel and stained by Ethidium bromide. The samples 
were extracted with MinElute Gel Extraction Kit (Qiagen, 
Germany), and these purified PCR products were cloned into 
pGEM-T plasmid vector (Promega, USA). Three clones were 
subsequently sequenced (Invitrogen, USA). The online tools, 
Open Reading Frame Finder (http:/ /www.ncbi.nlm.nih.gov/ 
gorf/gorf.html), Compute pI/Mw (http://web.expasy.org/ 
compute_pi/), and InterPro (http://www.ebi.ac.uk/interpro/), 
were used to analyze the nucleotide and amino acid sequences. 

Vector Construction and Plant Transformation 

The coding region of £mARF25wa$ introduced into pDONR221 
vector using Gateway BP clonase enzymemix (Invitrogen, USA) 
with the following primer pair: Gate-25-F and Gate-25-R (Table 
SI). To construct overexpression vector 35 S::£mARF2 5, the 
fragments were transferred from the pDONR221 vectors to 
pB2GW7 vector by Gateway LR recombination (Invitrogen, 
USA). This construct was transformed into Arabidopsis (Col-0) via 
the floral dip method using Agrobacterium tumefaciens strain 
GV3 101 [35] . Transgenic plants were selected on 1 /2 MS medium 
containing herbicides (Basta). The herbicide resistant seedlings were 
transferred to a mixture of soil and vermiculite (2:1) and generated 
homozygous lines by self-fertilization. Statistical analysis of the 
difference in traits between transgenic plants and wild-type plants 
was performed by using /-test. 

Results 

Comparison of Maize Leaf Area between Hybrids and 
their Parental Lines 

Previous studies revealed that maize hybrid exhibited heterosis 
in leaf size [24] . In this study, dynamic growth patterns of maize 
hybrid Mol7/B73 and its inbred parental lines were investigated 
in terms of the third leaf area and mid-parent heterosis (MPH) was 
calculated. The results exhibited that leaf areas of hybrid Mo 17/ 
B73 were significantly larger than its parental lines at 8, 10 and 12 
DAG. Remarkably, MPH (72.00%) at 8 DAG was much higher 
than that of 10 and 12 DAG (60.50% and 58.26%) [Figure SI]. 
Thus, the third leaf areas of eight hybrids and their parental lines 
at 8 DAG were further measured and compared. It can be seen 
that the leaf areas of all hybrids were significantly larger than those 
of their parental lines at 8 DAG [Figure 1A]. 

Identification of a Hybrid Down-regulated Gene 
ZmARF25 in Maize Seedling Leaves 

To gain insight into the relationship of ZjnARFlO and ZjnARF25 
with leaf size heterosis, expression patterns of these two genes in 
the third leaves of hybrids Mol7/B73 and B73/Mol7 and their 
parental lines. Considering of the time lag between differential 
gene expression and the observed heterosis of leaf size at 8 DAG, 
the third leaves at 5 DAG were selected for qRT-PCR analysis. As 
shown in Figure IB, the relative expression levels of £mARF25 in 
the basal region of Mol7/B73 and B73/Mol7 hybrids were 
significandy lower than the levels observed in parental lines, 
whereas no significant differences in the expression patterns of 
ZjnARFlQ were observed in these hybrids. Most notably, no 
obvious expression alterations were detected for ZjnARFlO and 
%mARF25 in the leaf tip regions of the hybrids and their parental 
lines [Figure 1C]. Furthermore, the expression patterns of 
ZjnARF25 between other hybrids and their parental lines were 
also investigated and we observed that ZjnARF25 was down- 
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Figure 4. Morphology of 35S::Zm/'l/?F25 and wild-type plants. Morphology of (A) leaf size, (B) Inflorescences, flowers and siliques and (C) plant 
height. Statistics of (D) fully expended fifth leaves area, (E) siliques and (F) plant height. Col represents the wild-type; #1-4 represent the four 
35S:ZmARF25 transgenic lines; data are represented as the mean ± SE (n = 10). f-test compared with the wild-type control: **refers to significant at 
p<0.01. 
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regulated in the leaf basal regions of Zong3/87-l, 87-l/Zong3, 
B73/87-1, 87-1/B73 and Zong3/B73, but not B73/Zong3 
[Figure ID]. Therefore, ZjnARF25 was selected for further 
investigation. 

To obtain the cDNA with complete open reading frame (ORF) 
corresponding to J(mARF25, cDNAs from B73 seedling leaves were 
amplified with ZmARF25-F/R primer pair and one unique PCR 
fragment of about 2500 bp was produced (data not shown). The 
PCR products were cloned into pGEM-T and three clones were 
sequenced. One 2552 bp fragment was obtained and no 
nucleotide polymorphisms were detected for all sequenced clones. 
Sequence analysis revealed that the PCR product contains a 
complete ORF, which encodes a polypeptide of 801 amino acids 
with a calculated molecular mass of 89.79 kTJ and isoelectric point 
of 6.09, which was consistent with our previous study [32]. 
Structure analysis exhibited that ZmARF25 has a typical ARF- 
type structure with a conserved DNA-binding domain (DBD) that 
consists of a plant-specific B3-type subdomain and an ARF 
subdomain required for efficient in vitro binding to the auxin 
response elements (AuxREs) in the N-terminus, a typical carboxy- 
terminal dimerization domain (CTD) with both III and IV 
subdomains and a middle region (MR) that might function as a 
repression domain (RD) since it is enriched in serine (S) [Figure 2] 
[36-39]. 

Functional Analysis of ZmARF25 in Arabidopsis 

To investigate the possible function of ^mARF25, overexpres- 
sion Arabidopsis transgenic plants were generated and herbicide- 



resistant To transformants were transferred into soil and grown in 
a greenhouse to obtain segregating T! progeny for genetic 
analysis. Of the 21 transformants, eight overexpression lines 
showed a 3:1 segregation pattern for herbicide resistance, among 
which four (35S::^jnARF25 #1-4) were selected for further 
analysis. As shown in Figure 3A, the insertion of the overexpres- 
sion vector was confirmed by genomic PCR. ZjnARF25 is a 
homologue of AtARF2 in Arabidopsis [32], the relative expression 
levels of exogenous gene £mARF25 and endogenous gene AtARF2 
in the transgenic and wild-type plants were analyzed by qRT- 
PCR. The results showed that different expression levels of 
ZjnARF25 were observed in these overexpression lines, but the 
transcription of AtARF2 in the transgenic lines was not affected 
[Figure 3B]. Subsequent analysis was conducted using these T 3 
homozygous transgenic plants. As shown in Figure 4, 
35S::J(/nARF25 #1-4 transgenic lines displayed markedly de- 
creased leaf areas, shorter siliques and plant height with smaller 
inflorescences and flowers in comparison with wild-type. Further 
analysis revealed that the size of the fifth leaves of 35S::£tnARF25 
plants was decreased by approximately 60~80% as compared to 
wild-type [Figure 4D] . In addition, the plant height and the silique 
length of transgenic plants were about 40-60% and 56% of that of 
wild-type, respectively [Figure 4E, F]. Collectively, these data 
indicated that overexpression of ZjnARF25 in Arabidopsis decreased 
the size of aerial organs by inhibiting organ growth and 
development. 

To determine whether it is cell number or cell size, or both that 
contribute to the changed leaf size in transgenic lines, we observed 
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and compared die size and number of cells in the fully expended 
fifth leaves between 35S::£mARF25 and wild-type plants. The 
results showed that there was no distinct difference in cell size 
between the £mARF25 overexpression and wild-type plants 
[Figure 5A, B], whereas cell number of 35S::£mARF25 plants 
was significantly lower when compared to wild-type plants 
[Figure 5C], indicating that the decrease in cell numbers, not 
cell size, was associated with the diminution of leaf size in 
35S::ZmARF25 plants. 

Expression Analysis of Organ Size-associated Genes 
between 3SS::ZmARF25 Transgenic Lines and Wild-type 

To elucidate the underlying regulation mechanism of ZjnARF25, 
expression patterns of organ size-associated genes AtANT, AtGIFl, 
AtGRF5 and AtEXPK) between 35S::Z>nARF25 and wild-type 
plants were analyzed by qRT-PCR. The results showed that the 
expression levels of AtANT, AtGIFl and AtGRFS were significantly 
decreased in 35S::J(jnARF25 plants, whereas no obvious alteration 
of AtEXPIO expression was observed [Figure 6]. In Arabidopsis, 
AtANT, AtGIFl and AtGRFS were reported to be involved in cell 
proliferation while AtEXPIO played important roles in cellular 
expansion [40-44]. Collectively, it could be concluded that 
overexpression of %mARF25 inhibited the transcription of cell 
proliferation associated genes, which is consistent with the 
decreased cell numbers observed in 35S::J(/nARF25 plants. 

Discussion 

ZmARF25 Repressed Organ Size by Decreasing Cell 
Numbers in Arabidopsis 

Auxin response factors (ARFs) are transcription factors that bind 
with specificity to TGTCTC-containing AuxREs found in 
promoters of primary/ early auxin response genes and mediate 
responses to the plant hormone auxin [36-39]. Genome-wide 
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Figure 5. Histological analysis of fully expanded fifth leaves of 

35S::ZmARF25 and wild-type plants. (A) Epidermal cells of the fully 
expanded fifth leaves of 35S::ZmARF25 and wild-type plants. The 
statistics of (B) epidermal cells size and (C) cell number per leaf. Fully 
expanded fifth leaves of wild-type transgenic plants were analyzed. 
Data are represented as the mean ± SE (n = 10). f-test compared with 
the wild-type control: **refers to significant at p<0.01. 
doi:1 0.1 371 /journal. pone.0094830.g005 




Figure 6. Expression patterns of organ size related genes in 

35S::ZmARF25 and wild-type plants. The aerial portions of 1 5-day- 
old seedlings were used for qRT-PCR analysis. ACTIN 2 was used as the 
endogenous control. Col represents the wild-type; #1-4 represents the 
four 35S::ZmARF25 transgenic lines. Data are represented as the mean 
± SE (n = 3). f-test compared with the wild-type control: *refers to 
significant at p<0.05, **refers to significant at p<0.01. 
doi:10.1371/joumal.pone.0094830.g006 

analysis revealed that ARF genes are represented by a large 
multigene family in plants [32,45-49]. Most ARFs contain a 
conserved N-terminal DBD, a variable MR and a typical CTD 
[36,37,45-49]. The DBD specifically binds to the conserved 
AuXRE, the CTD is modular with amino acid sequence related to 
domains III and IV in Aux/IAA proteins, making it function as a 
dimerization domain among the ARFs or with several Aux/IAA 
proteins. The structure of MR determines whether the ARF acts 
as an activator or repressor. Activation domain (AD) of ARFs is 
usually enriched in glutamine (Q), serine (S) and leucine (L), while 
repression domain (RD) is enriched in either S, L and proline (P); 
S, L and/or glycine (G) or S [32,39]. In this study, one hybrid 
down-regulated gene %rnARF25 that encoded a typical ARF-type 
protein was identified in maize. Most notably, ZmARF25 has a 
MR enriched in S, indicating that it might act as a repressor. Our 
previous study indicated that ZmARF25 was homologous to 
AtARF2 in Arabidopsis and OsARF4 in rice [32]. AtARF2 has been 
identified as a general repressor of cell division, one aspect of many 
biological processes regulated by auxin [29]. However, our 
knowledge about the functions and regulations of ZjnARF25 gene 
is quite limited. Here, we demonstrated that overexpression of 
J(jnARF25 in Arabidopsis resulted in significant size reduction of 
plant organs. Although the results of ^mARF25 overexpression in 
Arabidopsis may not fully reflect its precise functions in maize, it is 
reasonable to assume that ^mARF25 may be involved in maize 
growth and development process, which deserves for further 
investigation. 

At the cellular level, organ size is determined by two 
developmental forces,cell proliferation and cell expansion, which 
partially overlap in time and are assumed to be coordinated [50]. 
Nevertheless, larger organ tends to be composed of more cells 
rather than bigger cells in different species [5 1] . In present study, 
we noted that overexpression of 2jnARF25 in Arabidopsis signifi- 
cantly decreased the cell number, resulting in smaller organs in 
35S::£mARF25 transgenic plants when compared to wild-type. 
Recently, some regulatory genes that affect leaf development have 
been identified and characterized at the molecular level. For 
example, loss of AtANT function reduces the size of all lateral shoot 
organs by decreasing cell number [41]. AtGIF and AtGRF5 act 
together and are required for the development of appropriate leaf 
size and shape by promoting and/or maintaining of cell 
proliferation activity in leaf primordial [42,43]. On the other 
hand, expansins are thought to function in the control of plant cell 
growth, and leaf size was substantially reduced in antisense lines 
with suppressed AtEXPIO expression, whereas overexpression of 
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AtEXPIO resulted in plants with somewhat larger leaves [44]. 
Remarkably, overexpression of £mARF25 repressed the expression 
levels of AtANT, AtGRFS and AtGIFl, whereas no significant 
alteration in AtEXPIO expression was detected. Collectively, we 
speculate that ZjnARF25 regulates the process of cell proliferation 
by altering the expression of AtANT, AtGRF5 and AtGIFl in 
Arabidopsis. 

Expression Patterns of ZmARF25 in Maize Hybrids and 
their Parental Lines was Tissue and Genotype Dependent 

Leaf differentiation proceeds basipetally (from tip to base) in a 
highly regular and continuous manner [52-54]. Maize leaves have 
a basal growing region, which is composed of immature cells, 
enclosed by older leaves [19,33,55]. Dynamics of the third leaf of 
maize seedling transcriptome and proteome along the maize leaf 
developmental gradient have been discussed. The results showed 
that genes and proteins expressed at the highest levels near the 
base of the leaf involved in cell wall biosynthesis, DNA synthesis, 
cell cycle regulation and chromatin structure regulation, protein 
metabolism, potential signaling proteins, auxin and brassinosteroid 
biosynthesis and signaling, respiratory pathways and vesicle 
transport [19,56]. Interestingly, relative expression level of 
ZjnARF25 gene in the basal region of hybrids Mol7/B73 and 
B73/Mol7 was significantly lower than that of two parental lines. 
On the other hand, no significant alteration in expression was 
observed in the tip region, which is consistent with the observation 
that different regions within a leaf displayed different expression 
patterns [19,56]. Taken together, we proposed that leaf basal 
region was the appropriate tissue for further investigation of the 
molecular basis of leaf size heterosis. 

Recently, many studies have analyzed heterosis-associated gene 
expression in leaves of maize, rice and wheat [7-14,20], but there 
has been no obvious consensus about genes that are differentially 
expressed in hybrids, which might be the result of significant 
differences in species, different genotypes within species, distinct 
tissues, experimental designs and statistical procedures applied [7- 
14,20]. In this study, we noted that the ^}nARF25 gene was down- 
regulated in B73/Mol7, Mol7/B73, Zong3/87-l, 87-l/Zong3, 
B73/87-1, 87-1/B73 and Zong3/B73, but no significant expression 
difference was detected between B73/Zong3 and its parental lines, 
which may be attributed to the influence of maternal affect, because 
the mid-parent heterosis of hybrid B73/Zong3 was lower than that 
ofZong3/B73. 
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